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I report a study of the nonstationary one- dimensional Fokker-Planck solutions by means of the strictly isospectral method of 
supesymmetric quantum mechanics. The main conclusion is that this technique can lead to a space-dependent (modulational) 
damping of the spatial part of the nonstationary Fokker-Planck solutions, which I call strictly isospectral damping. At the 
same time, using an additive decomposition of the nonstationary solutions suggested by the strictly isospectral procedure and 
by an argument of Englefield [J. Stat. Phys. 52, 369 (1988)], they can be normalized and thus turned into physical solutions, 
i.e., Fokker-Planck probability densities. There might be applications to many physical processes during their transient period. 

PACS numbers: 05.40. +j, 11.30. Pb 



At present, supersymmetry is recognized as a basic 
symmetry of natural world being discussed with vari- 
ous degrees of sophistication in many research areas. We 
also know that in physics it happens quite often that 
simple mathematical procedures might be extremely ef- 
ficient and may lead to significant progress in clarifying 
our vision of the world. This is undoubtedly the case of 
Witten's supersymmetric quantum mechanics [jjj which, 
mathematically speaking, is the factorization of the one 
dimensional (ID) Schrodinger operator and relies on the 
Darboux transformations. The procedure involves Ric- 
cati equations for the so called superpotential and most 
people got used to work with the particular Riccati solu- 
tion. However, Mielnik Q realized that one can use the 
general Riccati solution and presented the ID harmonic 
oscillator in that perspective, hinting on the connection 
with the inverse scattering methods. A clear discussion 
was provided by Nieto || and later Sukumar M showed 
that the GcPfand-Levitan method can be interpreted as 
a sequence of two Darboux transformations. While the 
particular superpotential leads to the quasi-isospectrality 
of the partner Hamiltonians, the general one allows for 
the introduction of entire families of strictly isospectral 
of either "boson" or "fermion" partner systems. In this 
paper I shall call Darboux- Witten (DW) isospectrality 
this strict isospectrality, which is obtained by employing 
the general superpotential (other people call it the dou- 
ble Darboux method). In a recent paper, Sukhatme and 
collaborators S obtained bound states in the continuum 
by means of DW isospectrality and a couple of other ap- 
plications may be found in the literature In the fol- 
lowing, after quickly presenting the factorization of the 
ID Fokker-Planck (FP) equation, I shall tackle its DW 
(strict) isospectrality. Before proceeding, I recall that 
Darboux transformations (not DW ones) have been used 
by Englefield for a FP equation with a particular po- 
tential, while Hron and Razavy || studied some solvable 
models of the FP equation by means of the Gel'fand- 



Levitan method. The other two techniques similar to 
the DW isospectrality, namely, the Abraham-Moses pro- 
cedure H and the Pursey one |l(J can be applied to the 
FP equation without any difficulty. 

Bernstein and Brown |nj were the first people to pro- 
vide a simple discussion of the correspondence between 
the ID FP equation with an arbitrary potential and Wit- 
ten's supersymmetric quantum mechanics. They showed 
that the great advantage of the supersymmetric proce- 
dure for the FP problem is to replace bistable "bosonic" 
potentials with much simpler single well "fermionic" 
ones. Here, I shall use the Smoluchowski form of the 
ID FP equation with constant diffusion coefficient and 
potential drift 
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where /' = df/dx is the drift force up to a sign, i.e., / 
is the FP potential, and 7 is a free parameter whose ex- 
pression depends on the physical problem under consid- 
eration. Eq. (1) is not a Hamiltonian evolution of the so- 
lutions. Nevertheless it can be cast into the Schrodinger 
equation as follows. Any initial time-dependent solution 
V(x,t) will relax at asymptotic times to the stationary 
solution 



V (x) = exp[—yf(x)] 



(2) 



but to pass to probability densities P one should intro- 
duce a normalization constant A^ 1 as a factor in the 
right hand side of Eq. (2). The issue of normalization 
looks obvious but let me emphasize that to make it clear 
one should consider the FP probability current. This im- 
portant quantity occurs when one wants to turn the FP 
equation into a continuity one. It reads 



J(x,t) = -e-T^e^P]' . 



(3) 
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The stationarity is defined as the case(s) of constant J, 
and the constant is determined by the boundary condi- 
tions. The particular case of zero J means no periodic 
boundary conditions and leads to e 7 ^Po = C, i.e., the 
constant is C = 1/Nq- To get the common notion of 
probabilities (a set of positive numbers, each smaller than 
unity, summing up to unity, and mapped on an equivalent 
set of future, possible events belonging to the same class) 
one should take the normalization constant as follows 



N c = 



(4) 



where the integration limits are from zero to oo and from 
— oo to oo for the half-line and the full line cases, respec- 
tively. 

The evolution at intermediate times can be discussed 
conveniently by means of the celebrated ansatz 



(5) 



V{x,t) — ► Vo when t — * oo, that is tp(x,t) oc y/Po(x) at 
asymptotic times. The ansatz (5) turns the FP evolu- 
tion of V into a Schrodinger evolution for the amplitude 
function ip in imaginary time 



-57 = -H F pip 
at 



(0) 



where the FP Hamiltonian is a Hermitian and positive 
semidefinite operator. It is now easy to proceed with 
the factorization and the whole Witten scheme for the 
FP Hamiltonian. We write Hppi = A\A^^ with A\ = 

,7//2 = 



e 7//2 —yf/2 = _§_ 
dx dx 



i 7 /' and A 2 = e~^/ 2 £, 



hjf. I recall that the usual quantum-mechanical 
factorization is of the sort H qm = AA^ + e, where e is the 
so-called factorization energy, which is a parameter fixing 
the energy scale. A zero value of this parameter, as in 
the case of the FP equation, indicates that the factoring 
of the Hamiltonian is done with respect to the "ground 
state" wavefunction. One can easily see that the FP su- 
perpotential is Wo = In ^/ 7 Po and thus the mentioned 
FP factorization is indeed with respect to the FP ground 
state amplitude. Also, I remark that the FP superpoten- 
tial Wo = —\lf is proportional to the drift force, which 
is a quite well-known result. The superpartner Hamilto- 
nian will be Hpp2 = A 2 A\. The two FP Hamiltonian 
partners can be written as follows 



Hfp,\,2 — -j— £ — 51,2 



(7) 



with the Schrodinger potentials S\ % entering simple Ric- 
cati equations 



At this point, it is worthwhile to emphasize the con- 
nection between the FP solutions and the Schrodinger 
solutions. According to Eq. (5), the FP solutions can al- 
ways be decomposed into a product of a non-stationary 
Schrodinger solution and the ground state amplitude of 
the corresponding stationary Schrodinger equation. As 
remarked by Englefield 0, \/Po in Eq. (5) can be any 
positive solution of the stationary Schrodinger equation, 
i.e., does not have to be normalizable. This is a crucial 
point for our arguments below. 

Let us pass now to the strictly isospectral construc- 
tion. In principle, one can use any solution tp(x) of the 
stationary Schrodinger equation to perform this construc- 
tion but we shall use the ground state amplitude 
ipo(x) = \/V^. It is saying that the strictly isospectral 
Schrodinger potentials are given by 

SSZi^ <Po\ A) = S 1 (x) - 2[\n(J + A)]" 

where A is a real parameter, that mathematically is the 
arbitrary integration constant of the general Riccati so- 
lution, and 



J(x) 



(10) 



where c = for the half line case and c = — oo for the full 
line one, whereas the expression of ip gen is given below. 
In this construction the ipo solution is reintroduced in the 
spectrum by using the general superpotential solution of 
the Riccati equation which reads 
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The way to obtain Eq. (|ll|) is well-known [|2|]3|] and will 
not be repeated here. From it one can easily get Eq. (9). 
The new class of solutions is a one-parameter family (p gen 
differing from <po, by a quotient 



(x;X) 



J + X 



= (p (x)/M D (x) 
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In the FP context we need <p gen > 0, i.e., A S (—T(x), oo). 
The range of A will be further restricted by a normaliza- 
tion condition (see below). The main point now is that 
the denominator Md(x) looks like a space-dependent 
modulation. Its general behavior has been obtained in 
||p|. There is a strong damping effect of the integral J', 
both for the family of potentials as for the wavefunctions. 
The parameter A is just signaling the importance of the 
integral term. There is some modulation close to the ori- 
gin and for small A parameters. At higher A-values the 
damping nature gets extremely strong and Lp gen is going 
rapidly to zero. 

The general form of the FP amplitude solution using 
the general Riccati solution reads 
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Vge 



exp 



exp(--7/)/Mb (13) 



and thus one can sec the self- modulation of the solutions. 
These modulated solutions are stationary solutions of the 
following Schrodinger equation 



d 2 tp(x) 
dx 2 



(14) 



corresponding to FP equations which are strictly isospec- 
tral to the initial one. Exactly as in supersymmetric 
quantum mechanics, the strictly isospectral FP equa- 
tions can be understood as a one-parameter family of 
"bosonic" FP equations having the same "fermionic" 
partner equation. Its factoring operators are B\ = + 



gen 



and B 2 = -§z — W 9en . The amplitude ipi = B\tp 
satisfies the time-dependent Schrodinger equation Eq. (6) 
unless £2 substitutes S\. Using now the essential fact 
that l/ifgen ("ground state zero-mode") is the general 
solution of Eq. (6) again when S2 replaces Si (see sec- 
tion 7.1 in the review of Cooper, Khare and Sukhatme 
[0), and being positive can be employed in Eq. (5) in- 
stead of \fpQ whenever ipi substitutes ip, one gets the 
following FP solution 

Vi{x,t) = V -j n ■ (Brf) = A . (15) 

OX \ <Pgen J 

If one calculates the normalization integral for this solu- 
tion one gets 



Ni 



P 1 {x,t)dx={^ n )\ 
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For confining FP potentials, i.e., going to 00 at both 
asymptotic limits one has — > for any A and there- 
fore Ni is zero. On the other hand, ip gen can be nor- 
malized and of course adding any multiple of Pi will not 
change the normalization constant. In this way, one ar- 
rives at strictly isospectral FP probabilities of the type 



P 2 (a;,*) = kPi(x,t) + N 2 <p 



gen ) 



where the normalization constant N gen is J12| 



Nnen = 



¥> 2 qen( x ) dx = 



(17) 



(18) 



The constant k is arbitrary, restricted by the physical 
condition P 2 > 0. Moreover, N gen imposes further re- 
strictions on the A parameter. The correct range of A in 
any physical problem related to FP transients should be 

A e (-l(x),-i) u (0,00). 

The decomposition Eq. (17) of the transient FP prob- 
abilities appears to be quite general. The behavior of 
the transient probabilities is dictated by the parameter 
A, which is also essential in the normalization issue. On 



the other hand, since S^ — > Si for A — > 00, the decom- 
position Eq. (17) does not apply any more in that limit, 
and one should use the common decomposition Eq. (5). 
In other words, although the isospectral solutions are sta- 
tionary ones they are related to the transient period and 
disappear in the asymptotic stationary regime. The pa- 
rameter A looks like an effective "time" (damping) pa- 
rameter measuring their importance during the transient 
phases. 

Perhaps it is worthwhile to remark that in our nota- 
tions the initial FP equation reads 
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Substituting A2 by Ai gives a FP form of the heat equa- 
tion 
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since the operator in the brackets is d/dx. 

The strictly isospectral FP equations can be written 
by analogy as follows 



Wt V ^ = dx- 



B2 + -7/ 



P(x,t) 



(21) 



and the corresponding isospectral "heat" equation will 
be 



s 9 

di V ^ = d-x 
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P(x,t) 
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Finally, I recall that besides the DW strictly isospec- 
tral potential Eq. (9), there are four other distinct fami- 
lies of strictly isospectral Schrodinger potentials obtained 
from pair combinations of Darboux-Witten, Pursey, and 
Abraham-Moses procedures O]. The formulas corre- 
sponding to Eq. (9) for the four families are collected 
in the paper of Khare and Sukhatme All of them 

are of the Darboux type and one can use any of them 
in the FP problem. From the scattering point of view 
all these families have the same reflection and transmis- 
sion probabilities. However, they differ from each other 
in reflection and transmission amplitudes. 

In conclusion, I have shown in a simple and general 
manner that DW strict isospectrality may provide inter- 
esting insights in the case of FP probabilities. 
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